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Abstract 

Knowledge of surface forces is the key to understanding a large number of processes in fields 
ranging from physics to material science and biology. The most common method to study surfaces 
is dynamic atomic force microscopy (AFM). Dynamic AFM has been enormously successful in 
imaging surface topography, even to atomic resolution, but the force between the AFM tip and 
the surface remains unknown during imaging. Here, we present a new approach that combines 
high accuracy force measurements and high resolution scanning. The method, called amplitude- 
dependence force spectroscopy (ADFS) is based on the amplitude-dependence of the cantilever's 
response near resonance and allows for separate determination of both conservative and dissipative 
tip-surface interactions. We use ADFS to quantitatively study and map the nano-mechanical 
interaction between the AFM tip and heterogeneous polymer surfaces. ADFS is compatible with 
commercial atomic force microscopes and we anticipate its wide-spread use in taking AFM toward 
quantitative microscopy. 



INTRODUCTION 



High-quality-factor resonators are increasingly used in a wide variety of ultra-sensitive 
measurements of force [1], mass [2] and motion [3]. High Q means that the oscillator can 
respond with large changes in oscillation amplitude or phase, to very small changes in a per- 
turbing force. The dominate measurement paradigm for exploiting this enhanced sensitivity 
is based on driving the resonator with one pure tone near its resonance, while monitoring 
the response to the perturbation as a change of amplitude or phase, or even a frequency 
shift when a feedback loop is used to lock the phase. In the field of AFM this measurement 
paradigm has been used to create images of surfaces to atomic resolution [1] with very little 
back-action on soft and delicate material surfaces. However, in spite of its name, the atomic 
'force microscope' does not actually measure the force between the tip and the surface while 
imaging in these modalities. The determination of tip-surface force requires that one mon- 
itors the change in response as one slowly changes the height of the resonator above the 
surface [5T-TI0]. Such measurements have provided impressive three dimensional force- volume 
plots[ll, 12J but they are fundamentally limited by extremely long acquisition times which 
necessitate enormous effort to reduce drift artifacts[13], making three dimensional force maps 
impractical for most applications. 

In this article we demonstrate a new approach to dynamic interaction measurement 
based on actively modulating the amplitude of oscillation while measuring the amplitude- 
dependence of the resonators response. A variation of amplitude which is slow on the time 
scale of the oscillation period lends itself to an analysis of the perturbing force in terms 
of the amplitude dependence of two signals: one in-phase and one quadrature to the sinu- 
soidal motion. The amplitude dependence of these two force signals allows us to determine 
the tip-surface interaction directly as a function of the cantilever deflection at fixed probe 
height. We show how these two quadrature force signals provide the foundation for model- 
free determination of the conservative force versus distance curve, as well as an energy loss 
versus amplitude curve. These two force signals represent essentially all the information 
that is possible to obtain about the nonlinear perturbing force, from analysis of oscillator 
motion in the limited frequency band near the high Q resonance. After deriving the theory 
of ADFS, we test the concept both with numerical simulation and experimental study of 
heterogeneous polymer surfaces. The improved force resolution of ADFS allows us to crit- 

2 



ically analyze a commonly accepted model of the tip-surface interaction and it enables an 
alternative approach for characterizing the mechanical properties of polymer surfaces. 



RESULTS 

Dynamics and force modeling 

In air and in vacuum the continuum dynamics of a cantilever beam are usually reduced 
to a single mode harmonic oscillator described by a simple equation of motion [H] 

z + -^z + ul(z - h) = j^-F drivc cos (wdrivct) + -fc F tt> (z, z, {z(t)}) (1) 

where z is the position of the tip in the lab frame, h corresponds to the static probe height 
above the surface and ujq = \J k /m is the frequency of the first flexural cantilever resonance 
with quality factor Q, effective stiffness k and effective mass m. The two forces acting 
on the oscillator are the tip-surface force F ts and the external sinusoidal drive force of 
strength F drive . The tip-surface force F ts is in general a complicated function depending on 
the instantaneous tip position z(t), velocity z(t) and, for hysteretic forces, on the past tip 
trajectory {z(t)}. One typically decomposes the tip-surface interactions into an effective 
conservative part which only depends on the tip position z and a more complicated effective 
dissipative part [6], 

F ts (z, z, {z}) = F c {z) + F nc (z, z, {z}) (2) 
When driving on resonance, arrive = w , the tip motion z{t) is approximately sinusoidal [15J, 

z(t) « A cos(uj t + <f) + h, (3) 

with amplitude A, frequency ujq and phase lag ip with respect to the drive force, even under 
the influence of the significant tip-surface interactions. This sinusoidal tip motion is a result 
of the large quality factor (Q w 100 in air, Q w 10000 in vacuum) which guarantees that the 
energy in the cantilever oscillation is greater than the tip-surface interaction energy during 
one oscillation cycle. 

The observation that the tip follows a simple orbit in the z-z phase plane leads to a 
more compact description of the tip-surface force . Every point in the phase plane is part of 
only one sinusoidal orbit with a specific amplitude. Together with the assumption that the 
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tip-surface force does not depend on previous interaction cycles, every point in the phase 
plane can be mapped to a unique history of the tip motion which starts at the maximum 
distance from the sample surface during an oscillation cycle. Thus, the dependence of the 
tip-surface force on the past tip trajectory can be incorporated into the dependence on the 
tip position z and velocity z . Therefore, the interaction is completely described as a two 
dimensional function in phase space. Such a function is depicted in fig. [lji. The plotted 
model force is a van der Waals-Derjaguin-Muller-Toropov (DMT) force with an additional 
damping term that depends exponentially on the tip position z. This force model is given 
by 

6{ao+h - z) 10 (4) 
-f^ + f£*v^(;z-/i) 3/2 + 7oexp(^Ji z>h 

where H is the Hamaker constant, R is the tip radius, is the inter-molecular distance, 
7o is the damping constant, z^ is the damping decay length and E* is the effective stiffness 
of the tip-sample system. The assumed numerical values for these constants are shown in 
table 1. 

Tip motion with fixed amplitude covers only one orbit in phase space (solid black lines 
in fig. [lja,). Force measurement at only one amplitude does not reveal the full nature of 
the tip-surface interaction as illustrated in fig. ([l]>e) where conventional force- distance 
curves for different oscillation amplitudes are plotted. It is not possible to extract the 
exponential dependence of the damping from only one measurement at fixed amplitude. 
Several measurements with different amplitudes are required to cover a larger region in the 
phase plane and gain a deeper insight into the nature of the tip-surface force, especially 
regarding dissipative interactions which can have very different physical origin. 



Force reconstruction 



The tip-surface force could be considered as a time-dependent force acting on the oscillator 
and the force could then be reconstructed from higher harmonics of the tip motion [TBI HZ]- 
However, equation ^ implies that higher harmonics of the tip-motion are insignificant. 
Indeed, their amplitudes are typically below the detection noise floor and special force 
transducers[18j, high interaction forces[16j or highly damped environments [i~9] are required 
for their measurement. The high quality factor of the cantilevers resonance amplifies only the 
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Fourier components of the force near the resonance frequency u>o, whereas higher harmonics 
are sharply attenuated. F ts (uo) can be determined from the tip motion and knowledge of 
the calibrated transfer function G, and the drive force -F drive , 



F ts {uJo) = G (u )Z(u ) - F drive (5) 



where Z(cu) is the Fourier transform of the tip motion z{t). With equation (|3j) the complex 
Fourier component F ts {ojo) can also be expressed as two real- valued Fourier components that 
are in-phase and quadrature to the motion 

F- = \.[ F t ,(z(t),i(t)) (z(t) ~ h) dt (6) 

*«=?jfM*<«>.*<«)>(^^ (7) 

where T = 2tt/uq is the oscillation period. Fj and Fq can be measured with high signal 
to noise ratio, providing a foundation for highly accurate force reconstruction. Fj and Fq 
are usually measured for different probe heights h above the surface [5HTU| . However, their 
dependence on oscillation amplitude A can also be used for force reconstruction. Hence the 
name amplitude dependence force spectroscopy (ADFS). 

The component Fj is called the virial of the tip-motion and is only affected by conservative 
tip-surface interactions |20j. Since the conservative tip-surface force F c depends only the tip 
position z, we can rewrite equation ^ as (see Supplemantary Information) 

Equation ^ shows that Fj(A) is the Abel integral transform of the conservative part of the 
tip-surface interactional]. The Abel transform has a unique inverse transform with which 
we determine the conservative tip-surface force as 

z dz J ^ z 2 _ A 

Thus, equation ^ enables force reconstruction at fixed static tip-sample separation h, using 
only the oscillation amplitude dependence of Fj(A). 

For dissipative forces such a simple reconstruction is not possible since they also depend 
on the tip velocity z. However, the force component Fq(A) can be interpreted as the 



amplitude-dependent energy loss per oscillation cycle (see Supplemantary Information), 



E dis (A) = -2nAF Q (A) (10) 
The amplitude-dependence of the energy loss gives a signature of the dissipation mechanism 



and equation (10) allows for the quantitative exploration of the dissipative interaction at 



fixed probe height h, whereas previously a slow change of h was required [22J 



Probing the tip-surface interaction 

To rapidly acquire the data for ADFS we use a multi-frequency lockin measurement 
scheme [23] which combines the fast tip oscillation at the resonance frequency Uq of the 
cantilever with a slowly varying amplitude. The resulting beat-like tip motion is periodic 
with the beat frequency Atu (fig. |2p). The calibrated transfer function G(uS) of the cantilever 
is used to convert the measured tip motion spectrum z(u) to the force signals Fi and Fq. 
Due to the slow change of the oscillation amplitude, every oscillation cycle at the fast 
frequency ujq can be considered as an oscillation with fixed amplitude. Over a full beat 
period a large region in phase space is covered (fig. 2d) and the complete dependence 
of Fi(A) and Fq(A) on the oscillation amplitude is revealed. The measurement time for 
complete Fj(A) and Fq(A) is 2 ms in our experiments under ambient conditions. This high 
acquisition speed makes it possible to combine high resolution surface imaging with high 
accuracy force measurements in a single scan with normal scan rates. To minimize feedback 
artifacts while scanning we evaluate Fj(A) and Fq(A) only while the oscillation amplitude 
increases. Moreover, the high quality factor resonance guarantees a gradual change of the 
oscillation amplitude, so that we can determine the tip-surface interaction without amplitude 
or deflection jumps that are caused by instabilities which are inherent to quasi-static and 
other dynamic force spectroscopy methods[24j. 



Numerical simulations 

We begin with numerical simulations of the equation of motion Q with the tip-surface 
force defined by equation ^ and the numerical values for the force constants given in table 
1. We consider a standard cantilever in air with resonance frequency of /o = "o/W = 300 kHz, 
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quality factor of Q = 400 and stiffness of k = 40 N/m. We apply a drive force with frequency 
fo and sinusoidally modulate the drive strength at a frequency of Af = 500 Hz. Without any 
tip-surface interaction, this drive force yields a beat-like motion with a maximum amplitude 
of 25 nm. In the presence of the tip-surface force, we let the system reach a steady state 
and then collect data during two cycles of the amplitude modulation at a distance of 17 nm 
above the surface. 

From the simulated motion we extract Fi(A) and reconstruct the conservative tip-surface 
force F c using equation (|9]). The reconstructed curve is shown in figure ([3]), and it is in 
excellent agreement with the actual force used in the simulations. At the sharp transition 
point from the attractive to the repulsive regime small deviations become visible in the order 
of 0.7 nN. 



Experimental results 

In material science often the local mechanical properties of heterogeneous polymer ma- 
terials are of interest. One example of such a material is a blend of polystyrene (PS) and 
poly(methyl methacrylate) (PMMA) with surface topography shown fig. The differ- 
ent solubility of these immiscible polymers causes the formation of domains with different 
height. The higher structures in fig. [4£i are PMMA-rich, while the lower background matrix 
is PS-rich [25] . Scanning this sample we collect data at 1024 points on each of 256 scan 
lines in 17 minutes. Without any additional assumptions about the force or the motion, the 
conservative force curve is determined using equation Q. A typical conservative force curve 
at one image pixel is shown in fig. [4Jd. The curve shows an attractive force close to the 
surface with a sharp onset of repulsive force when the tip makes mechanical contact with 
the surface. 

The ability to rapidly measure a conservative force curve with such high accuracy at 
every image point enables a more detailed study of the local nano-mechanical interaction. 
The repulsive conservative force is generally described in terms of a power law of the form 

F(z) = E*(z-z min ) p + F min (11) 

where E* describes a local stiffness and z m ; n is the position of the force minimum F m i n . 
A fit of this force model to the reconstructed force at every image point reveals that the 
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interaction exponent has a mean value of 2.35 ± 0.3, on both the PS- and PMMA-rich 
domains (see Supplemantary Informationn). This observation suggests that the geometry 
of the contact is close to the ideal case of a cone indenting a flat surface, where an exponent 
of p = 2.0 is expected |26j. We attribute this behavior to the actual conic shape of the AFM 
tip. In contrast, traditional quantitative AFM methods assume the DMT model geometry 
of a perfect sphere indenting a flat surface where an exponent of p = 1.5 is expected|27j. 
Fixing for the whole surface, we generate a map of the stiffness factor as shown in fig. 2c. 
Typically dynamic AFM requires an order of magnitude difference in stiffness in order to 
show image contrast. The high force accuracy of our measurement allows for the resolution 
of a factor of three difference in stiffness. The PMMA-rich domains are stiffer than the 
PS-rich matrix, in agreement with previous measurements [28] and in agreement with the 
expected bulk modulus[29j. 



Local stiffness is a property typically extracted from the conservative part of the surface 
force. However, Dissipative interactions are also present, which was recognized early in the 
history of dynamic AFM [151 EQ] ED] . These dissipative interactions are far more complex as 
they depend not only on tip position z, but also on the velocity and the history of tip motion. 



To get a signature of the dissipation mechanism we use equation (10) to reconstruct the 
energy loss per oscillation cycle as a function of the oscillation amplitude A. To demonstrate 
this we scanned a heterogeneous blend of polystyrene and low density polyethylene (LDPE). 
Figure |5^ shows the surface topography of the blend where LDPE aggregates to form droplet- 
like structures. Figures [5)d and ^jp show typical conservative force and energy loss curves 
reconstructed with ADFS on each of the two domains. For both domains the dissipated 
energy in fig. [5b increases significantly in the repulsive region of the conservative tip-surface 
force. For LDPE we see that both the dissipation and the size of the interaction region 
are much larger than for PS. Moreover, on LDPE the conservative force stays in the net 
attractive regime. We can interpret these observations as resulting from greater positive 
and negative surface deformation as the tip approaches and retracts above LPDE. This 
interpretation is supported by the fact that the LDPE domains are not in the stiff glass 
phase since the glass transition temperature of LDPE is below room temperature [29] . 
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DISCUSSION 



We demonstrated accurate and high resolution measurements of the conservative force 
and the dissipated energy while scanning. A direct analysis of the measured tip motion 
allows us to gain insight into the exact nature of the tip-surface interaction and dissipation 
mechanisms without extensive modeling based on idealized assumptions about the geometry 
of the tip and the surface. Such insight is often necessary to understand the contrast 
mechanism in AFM images. ADFS is compatible with the current trend in AFM toward 
smaller, higher frequency probes [311 132] because it requires only a narrow detection band 
near the cantilever resonance. With future generations of high frequency force sensors we 
envision that ADFS will enable high-resolution force-volume measurement at video rate. 
The general idea of exploring amplitude-dependence can also be extended to other modes 
of AFM like frequency-modulated AFM, torsional shear force microscopy or frictional force 
microscopy. Beyond AFM, the ADFS concept may inspire the development of new sensing 
schemes, as the basics concepts presented here are applicable to many types of measurements 
which exploit the enhanced sensitivity of a high-quality-factor resonator. 

METHODS 

Numerical methods 

The numerical simulations have been performed using the cvode adaptive step-size 
integrator [33]. Special care was taken to treat the non-smooth transition between attractive 
and repulsive forces at z = in the force defined by equation (]4]). 

Experimental methods 

Polystyrene (M w = 280 kDa), poly(methyl methacrylate) (M w = 120 kDa) and toluene 
were obtained from Sigma-Aldrich and used as purchased. The polymers were dissolved in 
toluene at a concentration of 0.52 wt% and spin-cast on a silicon substrate. The sample was 
scanned with a Multimode AFM and a Nanoscope 4 controller (Bruker). The PS/LDPE 
sample (Bruker) was scanned with a Dimension 3100 AFM and a Nanoscope 3a controller 
(Bruker). Both AFM systems were used together with a signal access module and a multi- 
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frequency lockin analyzer (Intermodulation Products AB, IMP 2-32) for drive synthesis and 
signal analysis. All scans have been performed with Tap300Al-G probes (Budget Sensors) 
which were calibrated with a non-invasive thermal method [51] . 

The authors acknowledge financial support from the Swedish Research Council (VR) and 
the Swedish Government Agency for Innovation Systems (VINNOVA). 
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TABLE I. Force parameters. Numerical values of the force parameters used for the tip-surface 
force in the numerical simulations and for the force plotted in fig. 1 

Force parameter Value 
Hamaker constant H 3.28 • 10~ 17 J 

Tip radius R 10 nm 

Inter-molecular distance ao 2.7 nm 



Effective stiffness E* 



1.50 GPa 



Damping factor 7 



2.2 -10~ 7 ^ 



Damping decay length z. 



'7 



1.5 nm 
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FIG. 1. Oscillation amplitude-dependence of tip-surfaces forces. A model tip-surface force is shown 
as a function of tip position and velocity (a). The force is given by equation Q. The zoomed inset 
shows the region close to the surface where the force is non-zero. The black lines are tip orbits with 
different amplitudes. The force is only plotted up to the maximum amplitude. The corresponding 
approach-retract force curves (b-e) show a clear amplitude-dependence of the dissipated energy. 
For small amplitudes, the force curve is essentially given by the conservative part of the tip surface 
force. For large amplitude, the dissipative interaction dominates and the forces experienced on 
approach and retract change significantly. 
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FIG. 2. Tip motion in the time domain (a, c) and in phase space (b, d). For constant amplitude 
only one orbit is covered in phase space (b) and it is not possible to reveal the complex dependence 
of the force on z and z. With modulated amplitude the tip performs a beat-like motion in the time 
domain (c). The amplitude modulation is slow compared to the frequency of the tip oscillation. 
At the lower turning points where the tip interacts with the surface, the motion can be considered 
purely sinusoidal with constant amplitude. Over a full modulation period, a large area in phase 
space is covered and the tip experiences the full z- and i-dependence of the tip-surface force (d) . 
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FIG. 3. Force reconstruction from simulated data. The reconstructed conservative force from 
simulated data and the actual are in excellent agreement force used in the simulation. 
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FIG. 4. Conservative tip-surface force measurements on the surface of a PS/PMMA blend. The 
plateaus of the surface topography (a) have a height of circa 5 nm and are PMMA-rich, while 
the lower regions are PS-rich. The red point marks the position for which the reconstructed 
conservative force is displayed (b). The force has an attractive minimum where the tip makes 
mechanical contact with the surface. Beyond this contact point the tip indents the surface and 
experiences a rapidly growing repulsive force. From reconstructions at all points, a map of the 
stiffness factor is generated (c). Separate histograms of the stiffness distributions (d) for the PS 
(blue) and PMMA (red) domains. The PMMA-rich areas are a factor of three stiffer than the PS 
matrix. 
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FIG. 5. Conservative and dissipative tip-surface interactions on the surface of a PS/LDPE blend. 
The droplet-like structures of the surface topography (a) are LDPE-rich. The red (LDPE) and the 
yellow (PS) spots mark the positions for which the reconstructed conservative force (b) and energy 
dissipation curves (c) are displayed. 
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